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ABSTRACT 

In Part I of this series, we introduced the Spherical Collapse (SC) approximation 
in Lagrangian space as a way of estimating the cumulants £j of density fluctuations 
in cosmological Perturbation Theory (PT) . Within this approximation, the dynamics 
is decoupled from the statistics of the initial conditions, so we are able to present here 
the cumulants for generic Non-Gaussian initial conditions, which can be estimated to 
arbitrary order including the smoothing effects. The SC model turns out to recover the 
exact leading-order non-linear contributions up to terms involving non-local integrals 
of the J-point functions. We argue that for the hierarchical ratios Sj, these non-local 
terms are sub-dominant and tend to compensate each other. The resulting predictions 
show a non-trivial time evolution that can be used to discriminate between models 
of structure formation. We compare these analytic results to Non-Gaussian N-body 
simulations, which turn out to be in very good agreement up to scales where a < 1. 

Key words: cosmology: large-scale structure of Universe-cosmology: theory-galaxies: 
clustering- methods:analytical- methods: numerical. 



1 INTRODUCTION 

The large scale galaxy distribution can be used to study 
the origin and dynamics of cosmological fluctuations. The 
one-point statistical clustering of matter density fluctuations 
Sr, smoothed over scale R, is characterized here in terms of 
the reduced Jth-order moments £j(R) = ( <5r ) . Assuming 
that gravity is the dominant dynamical effect, the evolution 
of £,j(R,t) is completely fixed by the initial conditions (IC): 
£j(0,-R). The problem is that even in the simplest case of 
pressureless (collisionless) matter, we do not have the exact 
solution to the dynamical equations. Here we concentrate 
on the approximation given by perturbation theory (PT) 
which provides a framework to study small departures from 
the linear theory solution. The question we want to address 
is: what are the predictions of PT for a given set of (Non- 
Gaussian) IC? This is important because large scale galaxy 
surveys can then be used to discriminate between models of 
structure formation (e.g., Silk & Juszkiewicz 1991). 

We already have a good idea of how to answer this 
question for Gaussian initial conditions (GIC), both at lead- 
ing order, or tree-level (Peebles 1980, Fry 1984, Bernardeau 
1992), and also with higher-order (loop) corrections (Scocci- 
marro & Frieman 1996a, 1996b, Fosalba & Gaztanaga 1998). 
Comparison with observations and simulations made it nec- 
essary to develop PT also for the smoothed fields (e.g., Go- 



roff et al. 1986, Juszkiewicz et al. 1993, Bernardeau 1994a, 
1994b). 

The analytic calculations in the literature for Non- 
Gaussian initial conditions (NGIC) refer to the variance, 
£2, skewness S3, (see Fry & Scherrer 1994, FS94 hereafter) 
and the kurtosis, 5*4, (see Chodorowsky & Bouchet 1996, 
CB96 hereafter). These studies are important from the the- 
oretical point of view, but of little help in practice because 
they only apply to the unsmoothed fields. Numerical work by 
Weinberg & Cole (1992) provided interesting insights on dif- 
ferent aspects of some particular Non-Gaussian models but 
did not address the general question of the initial conditions 
(other works include Moscardini et al. 1993). Gaztanaga and 
Mahonen (1996), studied the case of strongly NGIC in N- 
body simulations as compared to GIC, showing important 
differences that could be used to discriminate models with 
current observations. 

Here we follow the path set up in Paper I (Fosalba 
& Gaztanaga 1998) and study the case of a generic Non- 
Gaussian model for the IC, under the spherical collapse (SC) 
approximation, which will allow us to present predictions for 
the smoothed fields. In §2 we set the problem of NGIC and 
recall the SC model. In §3 we present the predictions for a 
generic case of NGIC and compare them to N-body simu- 
lations and some previous calculations. We end up with a 
discussion and the final conclusions in 54. 
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2 NON-GAUSSIAN INITIAL CONDITIONS 

For a given clustering in the initial conditions, we would like 
to derive the final clustering using the exact dynamics that 
rule the evolution of the underlying field. To characterize 
the clustering, we concentrate on the study the statistical 
properties of one-point Jth-order moments rrij =< S J > of 
smoothed density fluctuations 8. In particular we will use a 
combination of moments, the so call connected moments or 
reduced cumulants £j, which of a given J carry statistical 
information independent of the lower order moments, and 
are formally denoted by a bracket with subscript c. These 
cumulants can be easily obtained by using the generating 
function method: 



d J 



h = (?)=h*('") 



(i) 



A dimensional scaling of the higher-order moments in terms 
of the second order one £ 2 = a 2 is given by the following 
ratios: 



Bj 



a J -J/2 ■ 



(2) 



For Gaussian initial conditions (GIC) in PT, it is more useful 
to introduce the hierarchical coefficients, 



sj = 4^= b ' &] 1 " 



J/2 



£ 



(3) 



as PT predicts them to be time-independent quantities on 
large scales, for GIC. These amplitudes are also called nor- 
malized one-point cumulants or reduced cumulants. We shall 
also use skewness, for 53 and kurtosis, for S4. 

As in Paper I, we start by expanding the density con- 
trast 5 assuming that it is small, as is usually done in the 
context of perturbation theory (PT) in Euler space, 



S(x,t) = 5i(X,t) + 5 2 (X,t) + 5 3 (X,t). 



(4) 



were we consider S n « <J n -i. The first term 61 = 81 is the 
solution to the linearized field equations. The second term 
82 ~ 8? is the second-order solution, obtained by using the 
linear solution in the source terms, and so on. 

In Fourier space, we write, 8(k) = ^2°° = . 1 8n(K), being 
6 n (k) the n-th order perturbative contribution. The latter 
is expressed as an n- dimensional integral over the kernels, 
F n , that encode the nonlinear coupling of modes due to the 
gravitational evolution, 

S n (k) = d 3 q 1 ...d a q 1 S D (k-q 1 ...-q n ) x 

x F n (q 1 ...qJ8 1 (q 1 )..Mq n ) (5) 

where So is the Dirac function and the kernels F n are given 
by symmetric homogeneous functions of q n with degree 
zero, that is, some geometrical average (see Fry 1984, Go- 
roff et al. 1986, Jain & Bertschinger 1994, Scoccimarro & 
Frieman 1996a). 



2.1 Linear Theory 

If we only consider the first term in the PT series, Eq.W, 
and the growing mode, the cumulants of the evolved field 
will just be: 



where D is the linear growth factor and ( 8q ) correspond 
to the cumulants of the IC. Consistently, the hierarchical 
ratios (see Eq.H) will scale as: 



Sj = Sj{0)/D J 



(7) 



were Sj(0) are the initial ratios. Note that this implies that 
the linear growth erases the initial hierarchical ratios, so that 
Sj — * 0, as time evolves (and D — > 00). 

In terms of the dimensional scaling, see Eq.pl, we have, 
Bj — Bj(0), so that the linear growth preserves the initial 
values. Note that if we want to do a meaningful calculation 
of these ratios or the cumulants, in general, we might need 
to consider more terms in the perturbation series, Eq.|4j, 
depending on the statistical properties of the IC, e.g., how 
they scale with the initial variance, which is typically the 
smallness parameter in the expansion of the cumulants. 

For GIC both Bj(Q) = Sj(0) = 0, and we have to 
consider higher-order terms in the perturbation series to be 
able to make a non- vanishing prediction. 



2.2 Gaussian Tree-Level 

The computation of the cumulants in PT dates back to Pee- 
bles (1980) work where the leading order contribution to 
the skewness was obtained making use of the second-order 
PT, e.g., F2. Fry 1984 extended Peebles analysis by making 
the connection between tree diagrams (or tree-graphs) and 
the perturbative contributions to leading order in the GIC 
case. With the help of this formalism he was able to obtain 
the leading order contributions for the three- and four-point 
functions making use of the 2nd and 3rd order kernels, F2 
and Fz, in PT. Furthermore, Fry found that, in general, the 
lowest order (tree-level) connected part that contributes to 
£j are of order 2(J — 1) in <5i. Note that this involves the 
cancellation of J — 2 contribution to the moment of order J, 
m.j. This is a property of the GIC for which all \8{\ van- 
ish for J > 2. For GIC the leading order contribution is only 
given by tree-graphs and is therefore called the 'tree-level' 
(see also Paper I). 

Bernardeau (1992) found the generating function of the 
one-point cumulants to leading order for GIC. Here (and 
in Paper I) we present a simpler derivation, inspired in 
Bernardeau's work, that can be extended to higher-order 
(loop) corrections and NGIC. 



2.3 The Monopole Approximation & The 
Spherical Collapse Model 

Given the kernels F n in Eq. [p| , we define the monopole con- 
tribution to F n as the spherically symmetric (angle) average: 



Ft =< F n >= c n /n\. 



(8) 



The monopole approximation, 8 is the one resulting by sub 
stituting F n by < F n >. Under this approximation, we ob- 
viously have, from Eq.|5|: 

8n(k) = ^s l (k)*...*s l (k), 

n\ 

where * means convolution (in Fourier space), so that, in 
real space 8(k) = J^^Li ^n{k) becomes: 



0) 
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S(x) = Y, 5n(x) = /[*!(«)] = J2 5 [ ^ (a;)r 



(10) 



Thus, the monopole contribution to the cumulants in PT is 
given by a local-density transformation Eq.JlOJ, whose co- 
efficients c n are to be determined by the kernels F n , which 
are found by solving the perturbative equations under the 
relevant dynamics (in Fourier space). 

We can now easily estimate all the 1-point statistical 
properties in the monopole approximation to PT. This can 
be done by using the generat ing function method Eq.flfl, 
with the field S given by Eq.JlOJ. When the initial condi- 
tions are hierarchical, the resulting expressions can be found 
in Fry & Gaztanaga (1993), who consider a generic local- 
density transformation and find, to leading terms in 07: 



S- A = 
S 4 = 

s 5 = 
+ 



Si C + 3c 2 + O(af) 

Sf + 12c 2 S3 C + 4c 3 + V2c\ + 0(of ) 

si c + 20c 2 si c + isc 2 sl c2 + (30c 3 + 12002)53 



5c 4 + 60c 3 c 2 + 6OC2 + 0{aT 



(11) 



where Sj are the hierarchical amplitudes in linear theory 
(which are given by the IC). If the IC are not hierarchical one 
might have to consider more terms to estimate the leading 
order (see §H). These arguments are valid for any dynamics 
and also for both Euler or Lagrangian space, they apply to 
any leading order calculation where Eq.flq] is valid. 

To estimate the PT contribution to the cumulants (e.g., 
c n ) it is not necessary to calculate the structure of the ker- 
nels F n from the exact solution to the field equations of 8, 
as in Goroff et al. (1986). Given the equations for the evolu- 
tion of the field one can determine Eq.llOl and therefore c n 
by just requiring the solutions to be spherically symmetric. 
In Paper I, it was shown that for gravity, the spherically 
symmetric solution to the dynamical equations is given by 
the Spherical Collapse (SC), whose solution is well-known 
[e.g., see §4.1 in Paper I]. Thus, c„ = u n , without need 
of estimating the kernels F n or any integral. In particular, 
the v n coefficients are derived by taylor-expanding the para- 
metric solution to the SC model around Si = (see Paper I, 
§4.1). The connection between the Gaussian tree-level in PT 
and the SC model was already pointed out by Bernardeau 
(1992, 1994a, 1994b) although, there, the tree-level ampli- 
tudes were derived using the rather complicated formalism 
of the vertex generating function instead of the density con- 
trast itself. 

For GIC we also showed in Paper I how the monopole 
approximation gives the exact leading order, which is given 
by the tree-level solution. We next will argue that it also 
gives the exact contribution to the tree-graphs that appear 
for NGIC. 



( S 4 ) = ((Si + 5 2 + . . .f ) = ( St ) +4 ( 5l& 2 ) + .... (12) 

Now, 5 2 is given by the 2nd order kernel F 2 in Eq.H: 5 2 ~ 
F 2 Sf so that the term / SfS 2 \ above involves an integral 
over F 2 ( <5f V which is zero in the Gaussian case, but gives 
in general: 

(Sl) = (5l) c+ w(5!) c (Sl). (13) 

Consider first the second term of the above equation. The 
only connected part of the above term that contributes to 
the F 2 ( <5i ) integral in the kurtosis, is just given by a tree- 
graph, connected with no loops, of the type: 

/ d 3 qid 3 q 2 S D (ki q 1 - q 2 ) F 2 (q 1 ,q 2 ) 

(5 1 (q 1 )S 1 (k 2 )5 1 (h,))(5 1 (q 2 )5 1 (k 4 )), (14) 

here again, like in the Gaussian case [Eq.[23] in Paper I], the 
last two factors do not have any dependence on the angle, 
q i q 2 , and the contribution of the kernel F 2 is just given 
by the geometric average, e.g., the number c 2 /2. This illus- 
trates that the local transformation, Eg. \ 1 (]J, given by the 
monopole, also accounts for the exact contribution from the 
tree-graphs in the Non-Gaussian case, e.g. , the contribution 
12c 2 S^ c to Sa in the local expression ([ll]). It also accounts 
for the direct terms from linear theory [first terms to the 
right hand of equations dill)], which are zero in the Gaus- 
sian case, and are not regarded as tree- graphs. 

Nevertheless, for Non-Gaussian IC the tree-graphs do 
not necessarily include all the leading-order contributions, 
unlike the Gaussian case. There are terms like ( 81 ) in 

Eq.O], which also comes from (SfSz) and contribute to 
the kurtosis Eq.|12|. The contribution involves an integral 
of F 2 ( df \ . In the monopole approximation, Eq.jTol, this 

term is simply given by C2/2 ( Si ) . While in the exact 
calculation, its value could have contributions from higher 
multipoles, depending on the particular form of the initial 
conditions ( 5i ) , which is needed to do the integration. 
FS94 have estimated this multipole integration for several 
Non-Gaussian models for the initial kurtosis, £ 4 = / Sf ) . 
They typically find (see their Table 1) that the coefficients 
I {£,4]/ £,4, which corresponds to the relative contribution of 
the higher multipoles, are smaller than unity, thus typically 
smaller than the monopole contribution c n /n\ > 1 (see be- 
low Eq.p3J). This is not surprising as correlations typically 
decrease with distance, thus reducing the contribution from 
non-spherical geometries. Thus, in general, on expects the 
monopole to be the dominant contribution. In the case of 
hierarchical initial conditions, / Si \ is of order of, while 

< <5i > ( Si \ is of order of. Thus in this case the monopole 
(local) approximation exactly accounts for all the leading- 
order terms. 



2.4 Cumulants for a Non-Gaussian Tree-Level 

On estimating the PT predictions for cumulants £j, Fry 
(1984) described the connection between tree diagrams (or 
tree-graphs) and the perturbative contributions to leading 
order in the Gaussian case. These Gaussian tree-level terms 
will also contribute to the estimation £j in the case of NGIC, 
as they come from the products of linear terms, Si , induced 
by the two-point function. But there are additional terms in 
the Non-Gaussian case. For example, consider the kurtosis: 



2.5 Smoothing Effects 

In §4.4 of Paper I, it was shown how to relate the smoothed 
cumulants of the evolved distribution with the smoothed cu- 
mulants in the initial one. The arguments presented there 
were general and could be applied to both Gaussian and 
Non-Gaussian initial conditions. In the case of a power-law 
power spectrum P(k) ~ k n , the smoothed variance is also a 
power-law 07 ~ R 1 ' 2 , where 7 = — (n + 3). We then have: 
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m-fm + sy' 6 ] 



(15) 



up to a normalization factor given by Eq.[39] in Paper I. 
Note that this final result as well as the general expres- 
sion agrees with Bernardeau (1994a) arguments, based on 
the vertex generating function, but they do not limit us to 
Gaussian initial conditions or the leading-order term. Here 
again, the vertex generating function (/(— t), corresponds 
to cumulants in Euler space, while our local-density rela- 
tion f{Si), applies to Lagrangian space. To leading order, 
they both give identical results for Gaussian IC, but they 
yield different results in general for Non-Gaussian IC or for 
higher-order terms with Gaussian IC. 



3 PREDICTIONS FOR NON-GAUSSIAN 
INITIAL CONDITIONS 

Here we study the more general case of evolution in PT from 
Non-Gaussian initial conditions (NGIC), where £j(0) 7^ 0. 
If we use S n ~ S±, the first perturbative contributions are: 



ZAt) 



+ 



dUAo) + J/2 d 2 D j+ Uj+i(0) + - 



(16) 



where S" G 



Sj i0 (see notation in §3.1) are the Gaussian 
tree-level amplitudes, 



+ 7 



S? 



34 

7 
60712 

1323 



+ 



7 2 
7+37 



(17) 



and so on, for a power-law power spectrum and a top-hat 
window (see also e.g., §5 in Paper I.) These amplitudes are 
intrinsically gravitational as the GIC ones are zero for J > 2. 
In the case of NGIC, it is more difficult to present a per- 
turbation series because it depends on how £j(0) scales with 
£2(0). In a statistically homogeneous distribution, £j(0) — > 
in the limit £ 2 (0) — > 0, so that we can write £j(0) — > [£ 2 (0)] a , 
with a — a[J]. So we may consider the general case: 

h(0) = Bj [? 2 (0)] alJ] . (18) 

Let us consider the different possibilities for a = a[J] in the 
above series (bm. When a > J — 1 the initial conditions 
are forgotten, as the leading-order effect of the evolution is 
still the hierarchical term, 5j , which dominates the evolu- 
tion. In the more general case, the contribution to £j from 
Non-Gaussian initial correlations of order larger than J are 
suppressed by powers of the growth factor D and they only 
become important at late times. When J/2 < a < J — 1 
we have quasi-Gaussian but non-hierarchical initial condi- 
tions. The evolution in £j has a dominant non-hierarchical 
term that grows as D , while the hierarchical term grows 
as Z) 2 '- 7 " 1 ' and may not become significant until f 2 ~ 1. 
Note that, as pointed out by FS94, there is an additional 
Non-Gaussian term that grows as D J+1 and, for the skew- 
ness (J = 3), it contributes directly to 5*3 = £ 3 /£ 2 so that 
S3 7^ S3 at all times. When a = J/2, as in the dimensional 
scaling where £j(0) ~ £ 2 (0) J , both terms that grow as D 4 
depend on £ 2 and their total amplitude is S3 + 3/2c?2 . If 
a < J/2 there are strongly NGIC that dominate the evolu- 
tion as far as £ 2 is small (i.e., on large scales). 



Here we will concentrate on the a — J/2 case, the tran- 
sition to the strongly NGIC. For the dimensional scaling 
which typically arises in topological defects models such as 
textures, we have, 



G = B j£,2 > 



(19) 



where the Bj amplitudes are constants and have been pre- 
dicted to be of order unity (see Turok and Spergel 1991, 
Gaztanaga and Mahonen 1996). 



3.1 The Spherical Collapse Model Results 

In the SC model, it is straightforward to work out the per- 
turbation expansion. We can proceed as in Eq.JllJ by just 
keeping track of the dominant terms, as we do not need to 
solve any additional equations. We represent the different 
orders in the expansion of the cumulants following the no- 
tation introduced in Paper I: 

O 2 = \ S2,i a\ = o-f + S 2 ,3 of + S 2 ,4 o~i -\ (20) 



Sj 



, Sj,i a\ 



+ Sj-i o, + Sjfi 



+ Sj,i a i 



(21) 



and we find for the first non-vanishing perturbative contri- 
butions: 



S2,3 
S2,4 

53,-1 
5*3,0 

53,1 



sf 



B, 



3 ~3^ "18 ^ 3 > + T~ 



+ 



_^1 
2 

4 0) = B 3 



Sf -2 



(#) a S? 

12 12 



? 3 ° 



Bt + 



B 4 



§1 

2 



Bi 



si 

6 



17 



" V 5 3 G ) 2 +g5f 



B 3 



+ 



2 5 3 G + 



(*a G r 



-4+3 58-" 



3 



Bl 



54,-2 
54,-1 



2 r 
Sf = B 4 



6 
(5 3 G ) 2 



B3B4. 



_ + si 

12 8 



B s 



4 5 3 G B3+[3-5 3 G ] B3B4 + 



^ 



B B 



54,o = 54 + 



3 + 7 S3 



(5 3 G ) 2 + ^5 4 C 



m 



+ 



+ 



-i-Stf 



(5a G ) 2 5 - G 
~6~ + 4 Si 



Bi 



6-4Sf + |(S s G ) J 



BtBi 
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sc 


Unsmoothed 




Smoothed 




NGIC 


7 = 


7 =-l 


7 = -2 


7 = -3 


Bj = l 


n = —3 


n= -2 


71= -1 


n = 


S2,3 


0.62 


0.29 


-0.05 


-0.38 


S2,4 


1.87 


0.74 


0.44 


0.98 


s 2,5 


3.36 


0.60 


-0.05 


-1.05 


53,0 


5.05 


4.21 


3.38 


2.55 


•%,! 


7.26 


3.91 


1.55 


0.19 


Ss,2 


23.53 


7.37 


1.18 


0.20 


•54.-1 


19.81 


16.14 


12.48 


8.81 


54,0 


85.88 


52.84 


28.31 


12.27 


54,1 


332.51 


128.51 


32.83 


2.70 



Table 1. Values of the higher-order perturbative contributions 
in the SC model from NGIC with Bj = 1 for the unsmoothed 
(n = —3) and smoothed (n = —2, —1, 0) density fields for a top- 
hat window and a power-law power spectrum. 



-3+fs 3 G + 



m 



si 
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Bi 
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3S 3 G 



1 - - <? G - 

1 6 ^ 3 



-,G\2 



B 3 Be, 



'iT.si 

18 6 



B B 



SV, 



y(0) 



B 5 



C 9 -"U C 

5 S 3 B 3 + — S 3 Ba 



B 3 B 5 



+ 



5 „g 1 
gS 3 -1 



B 6 



56,-4 

Se,-3 



-r(O) 



B 6 



20 Sf B 3 B A + 10 Sf B 5 



sl 

3 



B;i Be 



+ [S 3 G - 1] B 7 . 



(22) 

We have chosen to write the expressions as a function of Sj , 
the tree-level hierarchical coefficients for GIC, so that when 
Bj = 0, one immediately recovers the Gaussian results. This 
is possible because, as mentioned in previous sections, the 
Gaussian tree- level contains all the information of the SC 
model (see e.g., Appendix Al in Paper I). Note that this 
includes both the smoothed and the unsmoothed case, by 
just replacing the correct value of n or 7 P in 5 G . We stress 
that, unlike the case of GIC, odd powers in the linear vari- 
ance also contribute to the cumulants for generic NGIC as 
already pointed out in recent papers (see FS94, CB96). 

Note that in Eq.f22|, there are two types of contribu- 
tions. First, there are terms that arise as a result of the 
non-linear gravitational evolution alone (~ Sj), which are 
not coupled to the IC, Bj and thus which contribute in the 
same way as for GIC. Second, there are those terms coupled 



to the IC (~ Bj) which depend on the specific Non-Gaussian 
model for the IC. 

It is natural to introduce the non-linear dimensional 
ratios for the texture- like NGIC (see Eq.JJ), 



B' 



= 2J Bj ' i a i = Bj '° + Bj ' 1 °" ;1 + 



(23) 



Notice that the linear term (~ 0(<tj)) remains equal to its 
initial value (contrary to the leading order to the Sj ratios, 
see Eq.|H]), while the cr-corrections (i.e., the 0(e>"j) terms) 
only appear when non-linear gravitational evolution sets up. 
In particular, we have to keep J — 2 corrective terms in order 
to include the purely gravitational term (not coupled to the 
IC), ~ Sj . Note that one has to include this hierarchical 
term to get an accurate estimate of the Bj , as its contri- 
bution increases with decreasing scale (see Figs U] & H for 
Sj = Bj0 2 ~ J ). 
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Bi 
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+ 
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2- y S 3 G 
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Bi 



(si) 2 
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Us?) 2 -is, 



Bi 



^s^ 2 -(s?r 



i-lsi 



r* (^ r 



18 



B 3 B$ 



Be 



(24) 



Notice that for most of the cosmological models, quasi-linear 
scales (where PT applies) have an effective spectral index 
n e ff, associated to a power-law power spectrum, P(k) ~ 
k n " f f , within the range n e // £ [— 1, —2], for which 5* G ~ 3. 
This means that, for most of the models, the second term 
in B 3j i & B4 t i have a negligible contribution. 

Table |I| shows the results for the NGIC for Bj = 1 and 
different power indexes. These are to be compared with the 
GIC case in Table 2 of Paper I. 

Figs hi and ti show the departures from the tree-level 
amplitudes for the variance, skewness and kurtosis up to 
the 2nd and 3rd perturbative order respectively, as the lin- 
ear rms fluctuation grows. It is seen how the 2nd order con- 
tribution (first corrective term) for the variance exhibits a 
cancellation of non-linearities for n»-l which tends to dis- 
appear when the linear rms density fluctuation approaches 
unity, once the next perturbative order (second corrective 
term) is included, as shown in Fig 0. We interpret that this 
is due to the fact that, here, unlike the case for GIC, the 
initial conditions still dominate the 2nd perturbative con- 
tribution in the evolved variance and further (3rd order) 
contributions change this picture since gravitational evolu- 
tion takes over the initial conditions. Put it another way, 
it takes longer for the gravitational evolution to erase the 
trace of strongly NGIC. 

On the other hand, for the skewness and kurtosis we 
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SC for NGIC: 2nd. Order 
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Figure 1. Departures from the tree-level contributions as the lin- 
ear rras fluctuation grows, for the variance, skewness and kurtosis 
of the density field as predicted by the SC model up to the 2nd 
non- vanishing perturbative contribution (first corrective term) for 
different values of the spectral index: the dotted line shows the 
n = —3 (unsmoothed) case, and the short-dashed (n = —2), long- 
dashed (n = —1), dot short-dashed (n = 0), dot long-dashed 
(n = 1) depict the behavior for the smoothed density field. The 
solid line shows the tree-level values (or linear term for the vari- 
ance) as a reference. It is shown the case of Non-Gaussian initial 
conditions with Bj « 1. 



SC for NGIC: 3rd. Order 
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Figure 2. Same as Fig 111 when the 3rd perturbative contribution 
is taken into account. 



see a significant change from the Gaussian case that holds 
in the 2nd and 3rd order analysis: strong non-linear effects 
are found for any value of the spectral index, contrary to 
what happened for GIC where there was a characteristic in- 
dex n ~ 0, where non-linear effects vanished. This difference 
must be due to the strong Non-Gaussian character of the ini- 
tial conditions. In fact, by looking at the expressions for S3 
and 5*4 for dimensional NGIC (see e.g., Eq.p2| above) we see 
that, unlike the case of the variance, the higher-perturbative 
orders are completely coupled to the initial conditions (~ Bj 
terms), so that, even for the smallest scales, deviations from 
the PT predictions for Gaussian initial conditions (where 
gravitational evolution clearly dominates) must show up. 
This deviation from the Gaussian prediction on small scales 
is realized as a shift from the Gaussian values in the hierar- 
chical contributions (order zero in a, Sj,o). 



3.2 Comparison with the Gaussian Case 

We have explored the non-linear corrections to the 1-point 
cumulants of the density field, within the SC model, for 
Gaussian (in Paper I) and Non-Gaussian dimensional initial 
conditions separately. The aim of this section is to provide a 
direct way of comparing the predictions in both cases so as 
to get a further insight on how the choice of initial conditions 
changes the non-linear evolution. 

To simplify the comparison, we shall assume that the 
Bj coefficients (see e.g., Eq.pl), are all equal, which is 
roughly what is expected from analytic models for Topo- 
logical defects. We denote = Bj as the Non-Gaussian 
strength. While the results in Table 1 assume /3 = 1, we will 
now display the non-linear corrections as a function of 0. 
We shall concentrate on the coefficients of even corrections 
in powers of a (S24, Sjfi), as the odd corrections vanish for 
GIC. Note nevertheless that these latter could be dominant 
contributions when studying the overall non-linear effect for 
NGIC. Using Eq.[E2|, we find for top-hat smoothing and a 
power-law spectrum: 



S-2A 
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(25) 



The case of GIC is reproduced by just setting 3 — 0, while 
one could expect 3 ~ 1 for a defect model. Recall that 7 = 
-(n + 3). 

Figures 0-o show the deviations from the Gaussian per- 
turbative contributions in the Non-Gaussian model as a 
function of the Non-Gaussian strength 8, for several smooth- 
ing values, n. They show the robustness of the arguments 
drawn from the 8 — 1 case (see Table 1). The Non-Gaussian 
contributions are larger than the Gaussian one for 1 > 3 > 
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Figure 3. Departures from the Gaussian behavior in the coeffi- 
cient S2,4 of the <t 4 correction to the linear variance as a func- 
tion of the Non-Gaussian strength /3 = Bj. The values displayed 
correspond to n = —0.5 (dotted line), n = —1 (short-dashed), 
n = —1.5 (long-dashed), n = — 2 (dot short-dashed), n = —2.5 
(dot long-dashed), n = — 3 (short-dashed long-dashed). 



0. For positive NGIC (/? > 0), this trend is typically en- 
hanced when smoothing effects increase. 

However, for negative values of the Non-Gaussian 
strength j3, the behavior changes. It is found that the vari- 
ance falls below the Gaussian value, while the Non-Gaussian 
Sj coefficients may overestimate or underestimate the Gaus- 
sian prediction depending on the scales (effective spectral 
index) we are looking at. Typically it underestimates the 
Gaussian value on quasi-linear scales. 

The general behavior found for the variance (s2,4,) seems 
to be in qualitative agreement with the conclusions drawn 
by FS94 which pointed out that a larger (lower) variance 
than the Gaussian one is expected for models with positive 
(negative) initial skewness (see also Moscardini et al. 1993). 
But note that this trend changes when S3 < 3, since, then 
S2,3 < 0, which is the leading term in the variance (see 
Eq.S). This was not detected by FS94 because they did 
not include smoothing effects. Thus, for n > —8/7 ~ —1.14 
or 7 < —1.85, there is a change in this trend, with lower 
variance for more positive skewness. In observations and also 
in CDM models (see Gaztahaga & Baugh 1998) , the spectral 
index is n < — 1 in the weakly non-linear scales, going to 
n > — 1 on large scales. This will produce a characteristic 
change in the shape (see Figure H below). 

For 53,o, a critical index is observed at n = —8/7 fa 
— 1.14 (7 = —13/7) below which the Non-Gaussian hi- 
erarchical amplitude has a maximum as a function of f3 
which depends on the value of n (see inset in Fig ^). For 
n > —8/7 ~ —1.14 there is instead a minimum as a function 
of /3. These bounds could be of interest when interpreting 
the estimations of S3 from observations. Given a value of 



Figure 4. Same as Fig. 3 for the hierarchical contribution to the 
skewness, 53,0- The inset in the bottom shows a detail of the plot 
for a shorter range of the Non-Gaussian parameter /3. 
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Figure 5. Same as Fig. rafor the hierarchical contribution to the 

kurtosis, 54 o- 



the index 7, there is a maximum ( 7 > —13/7) or minimum 
( 7 < —13/7) in the non- linear coefficient at: 



A 



3 (20 + 77) 



8 (13 + 77) 
for which the ratio is: 



S3,0 
^3,0 



7(2192 + 1624 7 + 2457 2 ; 
32(442 + 3297 + 49 7 2 ) 



(26) 



(27) 



For 54,0 w e can also see some bounds (which depend on 
the spectral index) that could be of practical interest. 
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3.3 Comparison with N-body Results 

The above predictions are in good agreement with the N- 
body simulations depicted in Figs pl-H. The N-body re- 
sults (symbols with error-bars from 3 realizations) are from 
Gaztanaga and Mahonen (1996). A non- linear sigma model 
was used for the texture dynamics, which was stopped at 
as — 0.1. Density perturbations were mapped by 100 3 par- 
ticles to produce the initial conditions for a gravitational N- 
body simulation. The density fluctuations were then evolved 
by a P 3 M-code (Efstathiou & Eastwood 1981, Efstathiou et 
al. 1988) until they reach as — 1.0. 

We first fit a model for initial conditions £j(0) using 
the statistics of the initial density fluctuations in the sim- 
ulations obtained from a map (which assumes fi = 1) of 
the texture dynamics (non-linear sigma model) at as — 0.1 
(e.g., 10 expansion factors before now: os = 1.0). For the 
initial variance we use a numerical fit, which will be scaled 
with a — D to be the linear variance. We corrected this ini- 
tial values from the grid shot-noise, and the results are the 
average of 3 realizations, which explains why the SC predic- 
tions are not smoothed in the Figures. For the higher-order 
cumulants of the IC, a good match for R > 10/i -1 Mpc is 
given by the dimensional scaling £j = Bja J , with Bj ~ 0.5. 
This model for the initial conditions is shown as the upper 
dotted lines in the plots for the Sj amplitudes (only larger 
scales are shown for clarity, smaller scales are dominated by 
large shot- noise fluctuations). 

For n ~ — 1 (7 — —2), the prediction for Bj — 0.5 is, 

a 2 « a 2 - 0.02 af 
Bf L « 0.5 + 3.1 tr, 
B? L « 0.5 + 6.2 ai + 19.7 af 



(28) 

We remark that a PT prediction for B™ must keep at least 
J — 2 corrective terms to the linear prediction to include the 
(hierarchical) purely gravitational term, Sf (see Eq.|2J]). 
This means that to yield a meaningful prediction for the 
skewness and kurtosis we have to include one and two (non- 
linear) corrective terms, respectively, to the linear theory 
prediction. In practice, n (or, equivalently, 7) is a function 
of the scale R, so that the coefficients of the non-linear cor- 
rections above (e.g., S2,3) are a function of the local slope, 7. 
This slope is obtained numerically from the initial variance. 
For the variance, £2 = a 2 , we find in Figure M a good 
agreement between the the simulation output (for as = 1.0, 
open squares) and the SC prediction for the leading contri- 
bution to the non- linear variance: S2,3 (short-dashed line). 
The agreement extends up to the scales where the shot- 
noise dominates the correlations (beyond the point where 
a ~ 0.5). In particular, the negative contribution predicted 
for the variance in S2,3 is clearly in agreement with the be- 
havior measured in the simulations. The latter is related to 
the fact that, contrary to the Gaussian case (see Paper I, 
§5.2), the variance has a critical effective index, n = —8/7 
for which the non-linear contributions change sign (see Fig 
fil). This has nothing to do with the previrialization effect 
found for GIC (since the SC model is unable to account 
for this effect, see Paper I, §5.2) but rather with the fact 
that the variance is still dominated by the IC which push 
the variance towards lower values. The prediction for the 
next contribution to the variance S2,4 for the NGIC (long- 
dashed line) and the Gaussian IC (dotted line), shows that 




8 10 20 

R (Mpc IT 1 ) 



40 60 



Figure 6. The ratio of the non-linear to the linear variance a /a? 
in texture-like non-Gaussian models. Displayed are the predic- 
tions from the SC model including ^2,3 (short-dashed line) and 
S2,4 (long-dashed), which is dominated by the Gaussian contri- 
bution (shown as the dotted line). The N-body simulation for 
as = 1.0 (filled triangles) shows a non-linear variance below the 
linear one on all quasi-linear scales in remarkable agreement with 
the SC prediction for the first non-linear contribution S2,3- 



this term is dominated by the Gaussian contribution (that 
independent of Bj), with a behavior similar to the one found 
for GIC (see Paper I, §5.2), i.e., positive contributions to the 
variance on all scales. A clear departure of the variance in- 
cluding the first non-linear correction (s2,3), with respect to 
the prediction once the second correction is included (s2,4) 
suggests the point beyond which the perturbative approach 
should break down (see arrow in Fig ph. The second non- 
linear (32,4) correction to the variance, seems to be dom- 
inated by the tidal contributions missed in the SC model 
prediction, as displayed in Fig H. Note nevertheless the sim- 
ilarity of this effect to the one found in Figure 6 of Paper I. 
There, the first order (non-linear) correction to the variance 
in the APM simulations also gives a much better agreement 
than the second order correction. This might indicate that 2- 
loop calculations have stronger tidal contributions, so that 
the SC model is a poorer approximation for higher-order 
loops. 

On the other hand, in the plots concerning the hierar- 
chical amplitudes Sj, the lower dotted lines are the linear 
theory predictions, which only approach the non-linear re- 
sults at the larger scales, R ~ 40/i _1 Mpc (a < 0.1), where 
they do better than the GIC predictions, but are not as good 
as the NGIC predictions. 

There is an excellent agreement between the SC model 
predictions for Sj with those from the N-body simulations 
up to the point where the prediction including the Sj t i con- 
tribution significantly deviates from that up to the hierarchi- 
cal term, Sj,o (which includes the purely gravitational term, 
Sq). This means that for the skewness, 53 (kurtosis, S4), the 
SC result is reliable as long as the prediction including the 
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Figure 7. The hierarchical skewness, 53, for texture-like non- 
Gaussian models. The triangles show the initial conditions (erg = 
0.1), which are fitted well by the dimensional scaling 53 = -B3/0", 
shown as the upper dotted line. Squares show 53 for a later out- 
put: erg = 1.0. The SC predictions for the og = 1 output are 
shown as short-dashed (including the second order contribution) 
and long-dashed line (including the third order). The continuous 
line shows the corresponding tree-level PT prediction for Gaus- 
sian initial conditions. The lower dotted lines correspond to the 
linear theory prediction. 



2nd (3rd) and the 3rd (4th) perturbative contributions are 
in rough agreement. This point is reached approximately for 
a ~ 0.5, which is shown as an arrow in Figs, m & H. Notice 
that for NGIC, loop corrections enter at the same order (in 
a) as tree-level corrections. Note also the deviation from the 
Gaussian prediction on quasi-linear scales which shows up 
as a shift from the Gaussian values, as mentioned above and 
predicted in Eqfcd] (order zero in a, Sj,o)- 

Similar results are found for higher-order moments, 
which have also been used to check the initial correlations, 
Bj ~ 0.5 up to J = 6. 

3.4 Comparison with Exact PT 

The analytic PT calculations available in the literature for 
Non-Gaussian initial conditions refer only to the variance, a, 
the skewness S3 (see FS94) and the dimensional kurtosis, B4 
(see CB96) for the unsmoothed fields. We will now compare 
our results to those of the exact (non-local) calculation. 

3.4-1 The Variance in exact PT 

For the unsmoothed variance FS94 find: 



-" =<rf + 21G + jim +0 



W) 



(29) 



where £ 3 is the linear 3-point function: £ 3 = (6f) c = D 3 ^(0) 
and I [£3] stands for some multipole integral of the linear 
three-point function contributing to the one-loop order. In 
particular, we define, 



Figure 8. The hierarchical kurtosis, 54, for texture-like non- 
Gaussian models. Symbols and lines as in Figure pj. The new 
line in the plot (dot long-dashed) displays the SC prediction in- 
cluding the 4th perturbative contribution. In this case the initial 
conditions are fitted to the corresponding dimensional relation 
S 4 = B4/0- 2 . 
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with cos# = XXf/xx/. That integral term, I[£s], competes, 
at the same order, with the tree-level amplitude in the per- 
turbative expansion for dimensional IC. As a result, though 
the tree-level amplitudes remain (local) shearless, the tidal 
one-loop contribution cannot be exactly recovered by the SC 
model. In particular, from the SC model we found for the 
evolved variance, 



~ 1 
0-1 + 



si 

3 



1 ) B 3 a? + O (af) 



(31) 



where here Bj = t^j/af denote the dimensional amplitudes, 
and S3 = 34/7 (the tree- level skewness for GIC). Replac- 
ing the latter number in the expressions above, we see that 
our results within the SC model are able to recover all the 
terms given by the exact calculation up to the one-loop tidal- 
induced term dependent on the initial three-point function, 
e-9-, f-^Kl]- Note that SF94 typically find (see their Ta- 
ble 1) that the coefficients r\j = I[(,j]/£,j resulting from the 
multipole integration are smaller than unity, r\j < 0.2, thus 
typically smaller than the monopole contribution . As men- 
tioned above this is not totally surprising as correlations 
typically decrease with distance, thus reducing the contri- 
bution in non-spherical geometries. Furthermore, note that 
the SC results above provides with the smoothed general- 
ization of the FS94 results (up to tidal effects), e.g., : 



(S 2 ) = erf 



+ 



13 7\ a 
21 + 3 ) ?3 



+ 



(32) 



The general behavior found for the variance seems to be 
in qualitative agreement with the conclusions drawn by FS94 
which pointed out that a larger (lower) variance than the 
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Gaussian one is expected for models with positive (negative) 
initial skewness (see also Moscardini et al. 1993). But note 
that this trend changes when S3 < 3, i.e., when S2.3 < 
0, which is the leading term in the variance (see Eq. 123]). 
This was not detected by FS94 because they did not include 
smoothing effects. Thus, for n > —8/7 ~ —1.14 or 7 < 
— 1.85, there is a change in this trend, with lower variance 
for more positive skewness. In observations and also in CDM 
models (see Gaztanaga & Baugh 1998), the spectral index 
is n 5s — 1 hi the weakly non- linear scales, going to n > —1 
on large scales. This will produce a characteristic change in 
the shape (see Figure H). 

Within SC model the reason for this qualitative change 
(lower or higher variance for more positive skewness) lies in 
the combination of smoothing effects and the transformation 
from Lagrangian to Eulerian coordinates, which reduce non- 
linearities. Indeed, the negative term in S2,3 = (S3 /3— 1) B3, 
which makes the cancellation possible, comes from the La- 
grangian to Eulerian transformation. While smoothing ef- 
fects, can reduce the value of S3 and change the relative 
sign. 

Note the repeated appearance of a a critical index 
n ~ -8/7 « -1.14 (7 ~ -1.85) in the SC model predic- 
tions for the cumulants (see also Paper I) . At this index, the 
corrections to the variance are minimal for GIC and change 
sign for NGIC. Furthermore, the skewness coefficient S3,o 
also shows a critical behavior (see §3.2) for that particu- 
lar value. This critical index happens to correspond to the 
solution to the equation S3 = 3, which lies close to the ob- 
servational values of both S3 (see Gaztanaga 1994) and n 
itself (Gaztanaga & Baugh 1998) in the quasi-linear regime. 



3.4-2 Skewness & Kurtosis in exact PT 



+ 



B, 
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B4 B c 
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(36) 



From the SC approximation we find a similar expres- 
sion (see Eq.|24|), which after replacing S3 = 34/7, matches 
the PT result for the tree-level contribution while failing to 
reproduce the one-loop tidal terms. These involve some in- 
tegrals of the initial J-point functions, just as the case of 
the skewness commented above. Again, we expect this tidal 
contribution to generically yield a marginal net correction 
to the SC model prediction (r/j < 0.2). 

Note again that all the exact results mentioned above 
correspond only to the unsmoothed cumulants. In the present 
context, the local contribution can be easily extended for 
the smoothed case. The smooi/iedpredictions are implicitly 
given in equations ( |22[ ) and Table M and tend to yield smaller 
non-linearities, similar to what was found in Paper I for GIC. 



3.5 Predictions for the isocurvature CDM 
cosmogony 

As an interesting working example, consider the isocurva- 
ture CDM cosmogony presented recently by Peebles (1998). 
In the particular model presented Peebles used for the ini- 
tial conditions the following parameters: B3 = D3 ~ 2.5, 
B± = D 4 ~ 9.9 and 7 = -2e = -1.2. From Eqs.[||,|| 
we have that the leading order non-linear corrections to a , 
B$ L = S s a and Sf L 



S4CT , in terms of the linear rms 



fluctuation, 07, are: 
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3 

,NL 
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(7; + 0.55 Oi 
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9.9 + 98 a x + 560 a? 



(37) 



Similarly, for the skewness, FS94 give 
S3 



£3 , 34 , 10 Ci 26 (£ 3 ) 2 
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(33) 



with the tidal terms expressed as integrals ~ /[^J,/^] of 
the linear three- and four-point functions respectively. Re- 
calling the corresponding expression for the skewness in the 
SC model, 



S 3 = ^+S 3 G - 
ai 



3 



u # + (%- 



1 B 4 +...(34) 



we see again how the SC approximation exactly recovers all 
the local terms, but not the tidal or non-local integrals. Note 
that, in this case, the two tidal contributions to the reduced 
cumulants appear at the one-loop order as differences in the 
J-point functions, leading in most of the cases (see FS94, 
Table 1) to a marginal contribution to the SC prediction. 

Finally, we turn to the result derived by CB96, for 
the dimensional (non-linear) kurtosis, B^ L for an arbitrary 
Non-Gaussian density field. They find, 



B? L =B 4 + 54,1 ai + O (of) 



(35) 



being, 
Ba 1 
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assumed B 5 ~ 1.6(B 4 ) 3/2 ^ 50, and B 6 ~ 3.3(B 4 ) 2 ^ 300, 
following the x 2 distribution. These non-linear corrections 
are not very sensitive to the values assumed for B$ and Bq , 
as setting them to zero (a very conservative assumption) 
yields B^ L ps 9.9 + 26 01 + 250 af, which is anyway a very 
large correction to the linear theory prediction, 54,0 = 9.9, 
when at ~ 1. We stress that a PT prediction for B^ L must 
keep at least J — 2 corrective terms to the linear prediction 
to include the non-negligible purely gravitational term, S*j 
(see Eq.p4J). This can also be seen in Figs m & H where the 
SC prediction follows the N-body predictions up to 01 < 1 
only when the first and second corrective terms are included 
in the perturbative expansion of the skewness and kurtosis, 
respectively. 

Note that the Isocurvature CDM model has a much 
larger non-linear correction for t he variance than that of the 
texture model discussed in 



3.3 



(see Eq.gJ and also Fig g), 
because there, the spectral index, is lower (n ~ —1) than 
in the Isocurvature CDM cosmogony (n — 2e — 3 ~ —1.8, 
see Peebles 1998). According to the above results, the non- 
linear corrections to B^ L are very large (and positive) and 
linear theory is no longer a good approximation for 01 ~ 1. 
This is very similar to the situation presented in Figures VYtt 
for the texture model, where linear theory only works for 
R > 50 ft -1 Mpc, or a; < 0.1. What is more, the non-linear 
corrections to the Bj L are typically an order of magnitude 
larger than the linear values (see Eqs.[M.[h7j). 
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These non-linear corrections should be taken into ac- 
count when comparing the model with observations in the 
galaxy catalogues, as in Table 1 of Peebles (1998), where 
even at the larger scales 8 ~ 1 deg (R ~ 10/i _1 Mpc), one 
has <ji ~ 1. Thus, for the parameters shown in Eq.[B7|, the 
model seems to be incompatible with current observational 
constraints from galaxy catalogues. 



4 DISCUSSION AND CONCLUSIONS 

For Gaussian initial conditions (GIC), we found in Paper I, 
that the Spherical Collapse (SC) model gives the exact tree- 
level contribution to PT. It was shown that this contribution 
can be derived by means of a local-transformation of the 
IC, what is much simpler than the vertex generating func- 
tion formalism developed by Bernardeau (1994a, 1994b). It 
was also seen in Paper I that the SC model (in Lagrangian 
space) also gives an excellent agreement for the hierarchical 
amplitudes Sj in the loop corrections, as compared with the 
results derived by Scoccimarro & Frieman (1996a, 1996b) for 
the exact PT in the diagrammatic approach. 

We stress the importance of applying the SC approx- 
imation in Lagrangian space. There, the SC model is de- 
scribed by a transformation that only depends on the value 
of the linear field at the same point (what we call a local- 
density transformation). However, when going back to Euler 
space the density fluctuation (defined at a point) is normal- 
ized with a factor which is a function of the (non-linear) 
variance. Since the variance is a volume average of the two- 
point correlation function, this factor yields some non-local 
contribution to the cumulants (in Euler space). This non- 
local contribution is missed when introducing the SC model 
in Euler space directly, thus is not surprising that the pre- 
dictions for the cumulants in the SC approximation in Euler 
space are a poor estimation of those in exact PT, as the lat- 
ter are dominated by the non-local (tidal) effects (see Table 
A2 in Paper I) . 

For the predictions within the SC model (in Lagrangian 
space) for the hierarchical ratios, Sj, tidal effects partially 
cancel out as they seem to contribute roughly hierarchically 
to the cumulants (see Appendix A2 in Paper I) . This is also 
true for the SC model in Euler space, but the dominance of 
non-local effects in the cumulants there, yield significantly 
different Sj ratios to those in exact PT (see Appendix A2 
in Paper I). Smoothing effects do not alter substantially this 
interpretation (at least for a top-hat window). 

As the SC model can be easily extended to the 
smoothed fields, we were able compare the predictions for 
the higher-order moments from the SC model to those mea- 
sured in CDM and APM-like N-body simulations, and they 
turned out to be in very good agreement in all cases up to 
the scales where o~i w 1, supporting our view that the tidal 
effects only have a marginal contribution to the reduced cu- 
mulants. Furthermore, the break down of the shearless ap- 
proximation roughly coincides with the regime for which the 
perturbative approach itself breaks down, a ~ 0.5. That is, 
where the contribution of the second and the third pertur- 
bative order in the SC model are significantly different. 

For non-Gaussian initial conditions (NGIC) the SC re- 
covers all tree-graphs exactly, including all contributions 
given by the exact PT for the variance, a, skewness, 53, 



and kurtosis, Si, up to some non-local integrals, /[£j], in- 
volving J-point initial functions. These last integrals arise 
as a result of the coupling between the asymmetric initial 
conditions with the tidal forces. We argue that for the hier- 
archical ratios, Sj, these non-local terms are sub-dominant 
and tend to compensate each other (e.g., I[£j] < £./)• 

The measured higher-order moments in the N-body 
simulations with NGIC (with dimensional, texture-like, scal- 
ing) turned out to be in good agreement with predictions 
from the SC model up to the scales for which the pertur- 
bative series breaks down (e.g., see Figures g-|8j). As men- 
tioned in Paper I there is a critical index, n* £ [—1,-2], 
for which tidal effects vanish and the SC is a good approx- 
imation even for the variance. This might explain why we 
find such a good agreement for the variance as compared to 
the simulations (which have n ~ — 1 on weakly non- linear 
scales). This good agreement is found when using only the 
leading-order correction 0[<rf], and breaks down after in- 
cluding the next term, 0{af\ (see Figure h). This effect is 
similar to what we found for GIC when comparing the vari- 
ance between SC predictions and simulations (Figure 6 of 
Paper I). This could indicate that 2-loop calculations have 
stronger tidal contributions, so that the SC is a poorer ap- 
proximation for higher-order loops. 

The general behavior found for the variance seems to 
be in qualitative agreement with the conclusions drawn by 
FS94 which pointed out that a larger (lower) variance than 
the Gaussian one is expected for models with positive (neg- 
ative) initial skewness (see also Moscardini et al. 1993). 
But note that this trend changes when S3 < 3 i.e., when 
S2,3 < 0, which is the leading term in the variance (see 
Eq.p3]). This was not detected by FS94 because they did 
not include smoothing effects. Thus, for n > —8/7 ~ —1.14 
or 7 < —1.85, there is a change in this trend, with lower 
variance for more positive skewness (see §3.4). 

The SC model predictions for Sj show how the NGIC 
evolve slowly towards the (Gaussian) gravitational predic- 
tions but, even at as = 1, are still significantly larger. They 
show a characteristic minimum with a sharp increase in 
Sj with increasing scales, just as found by Gaztanaga & 
Mahonen (1996) in N-body simulations of NGIC. The lat- 
ter agreement strongly constrains the tidal contributions to 
the higher-order moments and gives further support to the 
domain of applicability of the SC model in PT. Thus, the re- 
sulting SC predictions for the NGIC show a non-trivial time 
evolution that can be used t o st rongly discriminate models 
of structure formation (see §3.5). 
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